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ABSTRACT 

We have obtained the first constraints on extended Lya emission at z ~ 1 in a sample of five radio 
galaxies. We detect Lya emission from four of the five galaxies. The Lya luminosities range from 
0.1 — 4 x 10 43 erg s _1 and are much smaller than those observed for halos around higher redshift radio 
galaxies. If the z w 1 radio galaxies are the descendents the z > 2 radio galaxies, then their Lya 
luminosities evolve strongly with redshift as ~ (1 + z) 5 . There do not appear to be strong correlations 
between other parameters, such as radio power, suggesting that this observed evolution is real and 
not an observational artifact or secondary correlation. We speculate that this evolution of luminous 
halos may be due to gas depletion (as gas cools, settles, and forms stars) accompanied by an overall 
rise in the mean gas temperature and a decrease in specifc star-formation rate in and around these 
massive galaxies. 

Subject headings: galaxies: evolution — galaxies: individual (3C 210, 3C 265, 3C 266, 3C 267, 3C 324) 



1. INTRODUCTION 

The bright end of the luminosity function is domi- 
nated today by giant elliptical and cD galaxies. The 
history of these massive galaxies traces the evolution 
of the highest peaks in the initial density perturbation 
spectrum, and thus potentially places strong constraints 
on paradigms of structure formation. The host galax- 
ies of high-redshift radio galaxies (z > 1; HzRGs) are 
likely the progenitors of these modern-day giants: the 
HzRG if-band Hubble diagram is well-fit by the 'pas- 
sive' evolution of a s tellar population with a high for- 
matio n redshift (e.g., iMcCarthvl 119931 : iDe Breuck et al.1 
l2002t ). and NICMOS continuum images of 0.8 < z < 1.8 
3CR radio galaxies show that the starlight distributions 
are rou nd and symmetr i c, mostly with i? 1/,4 -law ligh t 
profiles (jZirm et al.lll999tlZirm. Dickinson, k Devil2003f ). 
Furthermore, the inferred stellar masses derived from 
modelling of the spectral energy distributions (SEDs), 
including data from the Spitzer Space Telescope, define 
the upper enve lope of the mass func tion for all high red- 
shift galaxies (jSevmour et al.l [2007) . In addition, age 
dating studies of a few z ~ 1 — 1.5 HzRGs suggests that 
they formed the bulk of their stars at significantly higher 
redsh ifts, z ^> 2, and evolved fairly 'passively' there- 
after (Stockton, Kellogg, & Ridgway 19951: iDunlop et ail 
Il996t ISpinrad et al.1 119971 ) . Finally, radio galaxies at 
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z > 2 are general ly found in overdense, protocluster, 
environments (e.g..lVenemans et aT]|2007HKodama et al.l 
[20071 IZ irm et al.ll2008f ) again consistent with them mark- 
ing high density peaks in the dark matter field. 

Many of the z > 2 radio galaxies, including those 
within protocluster regions, a re surrounded by gi- 
ant (> 100 kpc) Lya halos dMcCarthv et all Il990l; 
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19971 IVillar-Marti'n et al.l fl999t iReuland et al 
20031) . These halos have line luminosities of sev- 
eral xlO 44 erg s _1 , and suggest the presence of huge 
gas reservoirs. Initially these halos were thought to be 
associated only with, and perhaps powered by, rare, 
luminous radio galaxies. However, the discovery by 
Steidel et al. (2000) of two giant Lya "blobs" associ- 
ated with a galaxy overdensity at z — 3.09 and only 
loosely associated with any galaxies with detectable UV 
continuum suggest that Lya halos may be relevant to 
the formation of the most massive galaxies in general. 
Subsequent deep Lya imaging and spectroscopy of 
this same protocluster field identified many lower 
luminosity, radio-quiet, spatially-extend ed Lya emitters 
(jMatsuda et alj|2004 ISaito et al.ll2008[ ). 

The origins of the Lya emitting gas remain ambiguous. 
Monolithic cooling of pristine, infalling gas would have 



an extremely short lifetime (r c 



10 4 yrs; e.g., Dey 



et al. 2005) and could not explain the prevalence of the 
halo phenomenon. However, if the ionized gas reservoirs 
were replenished via further infall and/or outflows from 
the central source the 'monolithic' scenario may still be 
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viable. If we assume that the halos are formed by a mix- 
ture of inflow and outflow, the morphology and evolution 
of the Lya emission should provide clues to the nature 
of that mix and by extension the galaxy formation pro- 
cess. Another possibility is that each Hydrogen atom is 
ionized more than once on average. 

Little is known about how these halos evolve with red- 
shift. As the host galaxy and radio source age spectrally 
and dynamically, the interaction between the halo and 
its embedded sources may alter in form. Between z ~ 4 
and z ~ 1, the host galaxy is largely assembled as a 
mature giant elliptical, and a hot young stellar popula- 
tion is no mo re than a minor contributor t o the total 
UV flux fe.g- IZirm. Dickinson, k. Devll2003fl . The cen- 
tral ionization source may shift from being dominated by 
young stars and starburst outflows to being dominated 
by the interaction with the radio source and photoion- 
ization by the central engine. Understanding the con- 
nection between young radio galaxies at z ~ 3 and the 
more mature objects at z ~ 1 requires similar data at 
all redshifts, including measurements of their emission 
line properties. Moreover, studies of the same emission 
lines at similar sensitivity can distinguish genuine evo- 
lutionary trends from the consequences of line choice or 
instrumental limitations. 

In this paper, we present Lya measurements of 5 ra- 
dio galaxies at 0.8 < z < 1.3 obtained using the slitless 
prism mode with STIS on HST. This article is organized 
as follows: we present the observations and data reduc- 
tion in the next section (§ ??), the results in Section ??, 
a discussion of these results in the context of theories of 
galaxy formation in Section ?? and finally a brief sum- 
mary of our conclusions in Section ??. Throughout this 
paper we assume a cosmological model with 57a = 0.7 
and Q m = 0.3 and H = 70 km s -1 Mpc -1 . 

2. OBSERVATIONS 

We used the Space Telescope Imaging Spectrograph 
(STIS) on-board the Hubble Space Telescope during 
20 December 2001 - 12 April 2003 to image 5 high- 
redshift radio galaxies lying in the range 0.8 < z < 1.3 
(GO#9166; PI A. Zirm ). The target g alaxies are selected 
from the 3CR catalog (|Bennettlll962f ) and all are known 
to have strong nebular line emission based on their op- 
tical spectra. At the redshift of the target galaxies, the 
Lya line falls in the near-ultraviolet (NUV) and thus re- 
quires space-based observation. Observations were made 
with the STIS NUV-MAMAs, ultraviolet sensitive array 
detectors, with the prism, to obtain very low dispersion 
(~ 30A pixel -1 ) spectra over the Lya line to produce, ef- 
fectively, a line image. The extremely dark UV sky back- 
ground and absence of readout noise in the STIS MAMAs 
enhance our sensitivity to low surface-brightness emis- 
sion. Our resultant Lya images have depth comparable 
to those of z ~ 3 objects in a few orbits, with much 
higher spatial resolution. 

The dispersion of the STIS prism mode at redshifted 
Lya at z = 1 is » 3000 km s -1 pixel -1 , sufficiently 
low that the prism-dispersed Lya will essentially form 
a monochromatic emission line image, with the galaxy 
continuum dispersed beneath it. In addition to the slit- 
less prism data we also took direct images using both 
the NUV-MAMA (which includes Lya) and the optical 
CCD to determine the zeropoint for the source spectra. 



The details of exposure times are listed in Table [T] 

We extracted the Lyman-a fluxes in a simple man- 
ner. Visual inspection of the prism and direct image 
data show that the 2D spectra are dominated by line 
emission. The structures observed in the direct im- 
age are reproduced when the prism is used rather than 
being smoothed over larger areas. 3C 265 is the ex- 
ception where both continuum and the CIV emission 
line are also visible in the 2D spectrum. However, to 
model the spectra, even in the case of 3C 265, we have 
taken a simple power-law continuum with both f\ oc A 
and oc A~ 2 (i.e., flat in f\ and /„ respectively). We 
chose these two extrema based on the empirical range 
of far-U V slopes observed for both star-form ing galax- 
ies (e .g.. iMeurer. Heckman. fc Calzettilll999[) and AGN 
(e.g., IVanden Berk et alJl2001t iTelfer et al.M2002j ). We 
used the tabulated dispersion relation for the prism 
and the direct NUV image of each galaxy to construct 
continuum-only models of the prism observations. We 
thereby assumed that the morphology of the direct im- 
age, which includes Lya , is representative of the spatial 
distribution of continuum flux. These models were scaled 
and subtracted, to minimize the residuals, to leave only 
the line emission. We note that the exact details of the 
continuum-subtraction should not significantly change 
the results. Lyman-a fluxes were then derived by per- 
forming photometry in apertures designed to include the 
residual flux on these continuum-subtracted images and 
masking the CIV flux in the case of 3C 265. The resulting 
Lya line fluxes luminosities are listed in Table [1] 

The raw prism data, the continuum subtracted data 
and the data with a fiducial halo model added are pre- 
sented in Figure [T] The fiducial model halo has L^ ya — 
1 x 10 44 erg s _1 and follows a King model profile with 
core radius = 10 kpc, similar to observati ons of Lya ha- 
los su rrounding z > 2 radio galaxies (e.g.. lReuland et al.l 
2003). The observed optical and near-infrared data are 
from HST WFPC2 and NICMOS respectively. 

3. RESULTS 

Of the five z ~ 1 radio galaxies targeted, four have 
significant or marginal detections of Lyman-a flux. 
3C 210 shows no line emission to a limit of I/L yQ ^ 10 42 
erg s -1 . Of the four with Lya , two are bright, 3Cs 
266 and 324 with L Lya = 3.6 and 2.0 x 10 43 erg s -1 
respectively, while 3Cs 267 and 265 are both a factor 
of ten less luminous (see Table [lj. The evolution of 
the Lya luminosity around powerful radio galaxies 
with redshift is shown in Figure O The stars are the 
galaxies from the current stud y while the circles are 
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1997 


; iReuland et al.l 


2003: IVillar-Martm et all 


2003 


; IZirm et al. 


2005 


: IVenemans et al.l 120071 



Villa r-Martm et al.l 12007ft and the squar es are radio- 



quiet extended Lya em itters ('blobs'; ISteidel et al.1 
I2000t iMatsuda et al.ll200l . The solid lines are 'by eye' 
power law evolution models with ^Lya oc (1 + z) a where 
a = 5.0 and 2.3. The radial extent of these z ~ 1 halos 
are about a factor of two smaller than the high-redshift 
halos. There is no similarly strong redshift correlation 
with radio power or stellar mass of the radio host for 
HzRGs in general (see Seymour et al. 2007). In any 
case, it is clear that the lOOkpc, 10 44 erg s _1 halos 
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TABLE 1 

Radio Galaxy Observations and Measurements 



Galaxy 


z 


RA 


Dec 


NUV Prism 


NUV Direct 


Line Flux 


Line Luminosity 






(J2000) 


(J2000) 


Exp. Time (s) 


Exp. Time (s) 


10 — 15 erg s _1 cm -2 


10 42 erg s" 1 


3C 210 


1.17 


08:58:09.9 


27:50:52 


2757 


260 


< 0.1 


< 1 


3C 265 


0.81 


11:45:29.0 


31:33:49 


1330 


360 


28.6 


2.3 


3C 266 


1.28 


11:45:43.4 


49:46:08 


5260 


360 


1.5 


36 


3C 267 


1.14 


11:49:56.5 


12:47:19 


2709 


240 


0.2 


1.0 


3C 324 


1.21 


15:49:48.9 


21:25:38 


3474 


360 


3.4 


20 
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Fig. 1. — 2D prism spectrograms of the 5 radio galaxies (panel a) 
and their continuum-subtracted (panel b; Lya only for all except 
3C 265, where CIV is also visible) data and how they would have 
appeared with a fiducial L-^ ya = 1 X 10 44 erg s — 1 halo (panel 
c). Panel d shows the observed optical image and panel e is the 
observed near-infrared image. North is up, East to the left and the 
major tickmarks are separated by 2 arcseconds. 

generally seen around powerful radio galaxies at z > 2 




Fig. 2. — Trend of Lya luminosity with redshift for the four 
radio galaxies presented here (stars), higher re dshift radio galaxies 
(circles) and radio-quiet Lyq halos (squar es; Stcidcl ct all 120001 : 
IMatsuda et~aT][2u04l; IDev et al.ll2005l: Igaito et al ll2008l) . The two 
solid lines are power-law evolution tracks of (1 + z) 5 and (1 + z) 2,3 
to guide the eye. 



are not present around their counterparts (which are 
not necessarily their descendents) at z ~ 1. 

At least two of the radio galaxies (3C 265 and 266) 
show Lya emission with significant angular extent (Fig- 
ure [3]). For these galaxies, we also plot the fiducial halo 
model in the same apertures for comparison. In addi- 
tion to the difference in overall luminosity, the lower red- 
shift halos also appear to have steeper surface-brightness 
profiles with smaller extent along with some constant 
surface-brightness emission associated with the UV con- 
tinuum. The sizes of the detectable line-emitting regions 
are only tens of kpcs as opposed to the higher-redshift ex- 
amples with radii close to 100 kpc. More specifically, for 
the well-detected case of 3C 265, the morphology of the 
extended line emission seems rather different from the 
halos seen at high redshift. For the other galaxies, our 
data are not deep enough to reveal much morphological 
information. For 3C 265, however, we can see that the 



in 



3C 267 




Radius (kpc) 

Fig. 3. — Radial profiles of the Lyo emission around the 4 radio 
galaxies with significant line flux. For 3C 265 the regions where 
CIV emission is seen were masked for these surface-brightness mea- 
surements. The solid lines in the 3C 265 and 3C 266 panels rep- 
resent the surface-brightness profiles for the fiducial model halo 
(L hyct = 10 44 erg s" 1 ). 
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Fig. 4. — Ratio of number of ionizing photons as inferred from 
the X-ray luminosity to the number inferred from the Ly-o lumi- 
nosity versus redshift for a sample of radio galaxies drawn from the 
literature (circles) and from the current work (stars). 

extended line emission (panel a in Fig. [I]) roughly follows 
the spatial distribution of the UV continuum (panel d) 
which is aligned with the double-lobed radio axis, per- 
haps indicating that this line emission is excited by the 
AGN. This in contrast to the higher redshift halos which 
generally show little or no UV continuum associated with 
the line emission. 

To investigate the energetics of the halos, we have 
used X-ray core measur e ments from the literature 
(ICrawford fe Fabianl IT9961: IHardcastle fe Worralll [I999L 
Fabian. Crawford. fe Iwasawal [20021: iPentericci et al. 
20021: iDerrv et alj|2003t iFabian et alj|2003t iScharf et al. 
2003t lOverzier et al.l l2005f ) and Lv-q luminosities to 
make two different estimates of the UV ionizing radiation 
in the radio galaxy environments. We assume a single 
ionizing photon results in a single Lya photon. For the 
X-ray data we assume a single power-law spectrum from 
the rest-frame X-ray to the UV to estimate the number 
of ionizing photons due to the AGN. In Figure 2] we 
show the ratio of the two ionizing luminosities versus 
redshift. While this plot is a bit sparse and consists of 
rather heterogeneous datasets, it seems that there may 
be a trend for the high redshift halos to require a power 
source beyond the X-ray luminosity while at z ~ 1 the 
X-rays may be sufficient to power the extended line 
emission. We discuss the implications of this result 
further in the following section. 

4. DISCUSSION 

We have presented slitless NUV spectroscopy of five 
z ~ 1 powerful radio galaxies. These spectra show 
for the first time the considerable redshift evolution of 
the extended Lya luminosity surrounding powerful radio 
galaxies. The morphology of the extended line emission 
also seems to change from a centrally-concentrated halo 
at high redshift to a more evenly distributed surface- 
brightness profile at lower redshift associated with the 
known (extended) UV continuum. The size of the ha- 
los are also about a factor of two (at minimum) smaller 
than the high-redshift examples. These changes together 



suggest that we may be seeing a change in the dominant 
process responsible for the line emission as a function of 
redshift. Furthermore, the lack of similar halos around 
comparable power radio galaxies, along with the discov- 
eries of radio-quiet Lya halos, implies that halos are not 
solely a feature of radio galaxies but may be associated 
more generally with galaxy and structure formation. 

By making observations of Lya emission over a range 
of redshifts we study the evolution of the gaseous en- 
vironments of massive galaxies where other methods, 
such as X-ray imaging, are impractical at higher red- 
shifts. The luminosity and structure of the line emit- 
ting gas provides clues to the dynamical and excitation 
states of the gas. We can use simulations to put these 
observations in context. Dijkstra, Haiman and Spaans 
(2006a, b) have quantified the relationship between differ- 
ent infall and cooling scenarios and Lya surface bright- 
ness using a set of simulations. Their range of initial 
conditions and assumptions span the gap between two 
extremes. First, that the cooling timescale is short and 
tracks the loss of gravitational potential energy as the gas 
falls into the center after being shock heated at th e virial 
radius (e.g.. lHaiman. Spaans. fe QuataertJ l2000h . Sec- 
ond, that the gas is accreted cold (i.e., no virial shock 
heating) and only subsequently its kinetic energy is con- 
verted to thermal energy resulting in Lya emission (e.g., 
iDekel fc Birnboimll2006h . These simulations confirm and 
quantify the intuitive view that the line emission is more 
centrally-concentrated for the case where heating and 
cooling only occur once the gas has collapsed to the 
center and that more extended emission is seen where 
clumps of gas are assumed to cool as they descend into 
the gravitational potential. It should be noted, however, 
that both scenarios produce extended halos of line emis- 
sion. 

In addition to either or both of these infall scenarios the 
gas may also be ionized by UV-luminous young and form- 
ing stars and by the central AGN. For 3C 265, the close 
correlation between the UV continuum emitting regions 
and the Lya emission suggests that the ionizing source is 
low-level star-formation, young stars and photoionization 
by AGN light in these areas. Our observations indicate 
that the halo phenomenon as seen at high redshift has 
changed significantly by z ~ 1. This inference is further 
supported by our analysis of the energetics of the radio 
galaxies over the same range of redshift. While at high 
redshift z > 2 the virial shocking of infalling gas may 
dominate the Lya emission it seems that by z ~ 1 the 
X-ray output of the central AGN may be sufficient to 
ionize the surrounding gas. 

Under any scenario, it is clear that the gaseous envi- 
ronments of massive galaxies are subject to several pos- 
sible changes during the more than three billion years 
between z ~ 2.5 and z ~ 1. Gas is consumed in star- 
formation and black hole accretion, it is enriched with 
metals and deposited back in the IGM via supernovae 
and AGN-driven outflows and the IGM itself is disrupted 
and shock-heated via galaxy interactions and mergers. 
Overall, however, the amount of cool gas surrounding 
galaxies decreases with cosmic time while elliptical galax- 
ies and galaxy clusters retain vast reservoirs of hot gas. If 
we associate the high redshift, luminous halos with the 
initial accretion and starburst, the faint halos at later 
times suggests that the gas has either been depleted via 
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star formation and accretion or that the cooling time 
has increased sufficiently to quench infall. Cooling radi- 
ation and subsequent star-formation dominates the halo 
luminosity at high redshift with some additional contri- 
bution associated with the AGN. As the gas cools and 
forms stars the bulk of the line luminosity disappears 
with only the residual, AGN and star-formation related 
emission remaining by z ~ 1. This residual emission is 
over a smaller area (radius < 50kpc) and is associated 
with the UV continuum emission. 

5. CONCLUSIONS 

We have used STIS NUV-MAMA slitless spectroscopy 
of five z ~ 1 powerful radio galaxies to determine their 
Lya properties and to specifically look for luminous, ex- 
tended halos of line emission as are seen around z > 2 
radio galaxies. While we find that at least two of the 
targets have extended emission line regions it is unclear 
from these data alone whether this Lya is directly re- 
lated to the high-z halos. The emission line morphologies 
are rather different and the luminosities are an order- 
of-magnitude lower than at higher redshift. The physi- 
cal extent of the halos are smaller by about a factor of 
two compared to the z > 2 examples. We have argued 
that based on the available X-ray data it is likely that 
the dominant mechanism for producing Lya photons has 
changed from one related to gas infall and cooling at high 
redshift to one dominated by photo-ionization by young 
stars and the AGN by z ~ 1. 

Based on these data it is clear that Lya stud- 
ies of massive galaxies over a range of redshifts will 
provide clues to the state of the gaseous environ- 



ment and the history of galaxy feedback in the lo- 
cal IGM. In the scena ri o of Dekel and Birnboim 
(IDekel fc Birnboiml l2006t iBirnboim. Dekel. fc Neisteinl 
120071: iDekel fc Birnbo im 2008) there is a transition from 
a cold accretion regime, in which Lya photons are pro- 
duced by cooling gas, to a 2-phase medium, where Lya 
may still be emitted by filamentary 'cold flows', and fi- 
nally to a single hot phase. These transitions occur natu- 
rally as a result of the gravitational build-up of structure 
and the consequent formation of a virial shock. It may 
therefore be possible to follow this progression via Lya 
imaging surveys, particularly using high spatial resolu- 
tion data as we have presented here to study the morphol- 
ogy of the line emission. The state of the Lya emitting 
gas would provide direct clues to the state of the galaxy 
formation process and the thermodynamic properties of 
the gas. This " Lya calorimeter" could be combined with 
data from the next generation X-ray and mm/radio tele- 
scopes to generate a complete census of the gas in and 
around massive galaxies from z ~ 4 to z = 0. 
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